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Abstract 

I construct a minimal formalism to describe transport of non separable correlations 
arising when a superconductor is in contact with several electrodes (being ferromagnetic 
or normal metal). Transport theory is expressed in terms of effective single site Green's 
functions. Part of the circuit is decimated exactly by renormalizing only one physical 
parameter (the density of states). I show that the physical current is obtained by acting 
on the infinite series of Feynman diagrams with an operator counting the charge of a 
given diagram. I use this method to propose a new Aharonov-Bohm experiment intended 
to detect non separable correlations. In this experiment, spin-a electrons are forced to 
travel around an Aharonov-Bohm loop in the presence of an applied magnetic flux, while 
spin-(— a) electrons are not directly coupled to the magnetic flux. It is shown that the 
spin-(— a) current oscillates as a function of the flux to which spin-o" electrons are coupled. 
It is predicted that the effect exists both with high and low transparency contacts, and 
that the electrodes can be normal metal in which phase coherence can propagate over 
large distances. Three other aspects of the problem are investigated: (i) The dependence 
of the superconducting gap on spin polarization of ferromagnetic electrodes connected 
to the superconductor; (ii) The detection of linear superposition of correlated pairs of 
electrons; (iii) The proximity effect version of cryptomagnetism. 
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1 Introduction 



There are two important challenges in phase coherent transport phenomena. One is to fabricate 
quantum bits for quantum computing [I], g, g, ||, || || . The other is to make an Einstein, Podolsky, 
Rosen (EPR) experiment, that would be able to probe Bell inequalities with electrons 0]. Bell in- 
equalities have already been tested experimentally with photons [||, being massless particles. EPR 
experiments constitute one of the ultimate tests of quantum mechanics, which is a non local theory 
without hidden variables (at least with photons). At the present stage, there exists no clear proposal 
of what should be an EPR experiment with electrons. It is not even known whether it is possible 
to propose a realistic situation in which a transport experiment could detect the presence/absence 
of hidden variables with electrons. However, electrons carry an electric charge while photons do not. 
One can take advantage of this, and propose experiments with electrons that would not be possible 
with photons. 

It is already well established that the basic condition required to design an EPR experiment with 
electrons is to have a source of correlated pairs of electrons || Several dc transport experiments 



have been proposed recently |11, 12|, and a new proposal is made here. From a technical point of 
view, I concentrate on the simplest non perturbative formulation of transport theory, which relies 
on the use of effective single site Green's functions. By non perturbative, I mean that transport 
theory is exact to all orders in the tunnel amplitude. As a consequence, I can solve the physics of 
high transparency contacts. The interesting aspect of the formalism is that part of the circuit can 
be decimated exactly. Only one physical parameter is renormalized (the density of states). I show 
how to cure a more technical problem, related to the fact that all possible Feynman diagrams are 
present in the theory and have exactly the same value. The solution is to act on the infinite series of 
Feynman diagrams with an operator counting the charge of a given diagram. 

I apply the method to discuss four relevant physical questions: 

(i) What is the thermodynamics of a superconductor in the presence of correlated pairs of elec- 
trons ? I show that the strength of superconductivity depends on whether there exists or 
not non separable correlations in the vicinity of the superconductor. This may have direct 
implications for experiments. 

(ii) It is possible to fabricate linear superpositions of correlated pairs of electrons. How should 
one interpret transport of such states ? I show that transport theory is related to quantum 
measurement via projections of the linear superposition onto Cooper pair objects. 

(iii) Is there any way to manipulate one electron making the correlated pair and look at the response 
of the other electron ? I show that the answer is yes and propose a new type of Aharonov-Bohm 
effect that can be used to probe non separable correlations. Generically, there is no way to 
manipulate one electron without avoiding a response of the other electron. 

(iv) Is there a specific physics associated to Cooper pair penetration in ferromagnets in the presence 
of domain walls ? A simple model is solved, which shows that superconducting correlations 
can propagate along ferromagnetic domain walls. This may apply to recent experiments in fer- 
romagnet - superconductor heterostructures [13, 14, 15], and constitutes the proximity version 



of cryptomagnetism. 

Regarding the third question, it is worth recalling that in ordinary mesoscopic physics Aharonov- 



Bohm experiments [16, 17], electrons are forced to travel through a metallic ring in the presence of 
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a magnetic flux </> inside the hole of the ring. The current oscillates periodically as a function of <fi, 
with a period equal to the flux quantum cpQ = h/e |l8|]. This is because the system behaves like an 
interferometer in which the phase of the electronic wave function is related to the circulation of the 
vector potential. These oscillations occur in a variety of situations. For instance, the dual situation 
in which a magnetic flux is forced to travel around a magnetic charge is known as the Aharonov- 



Casher effect [19]. It is also very often that excitations of correlated states of matter do not carry a 
charge e. For instance, Cooper pairs carry a charge 2e. The associated flux quantum is h/(2e) [p0| . 
Now when ferromagnetic or normal metal electrodes are connected to a superconductor, it is possible 
to fabricate correlated pairs of electrons of the type c^c^JO) in which the spin- up and spin-down 
electrons making the Cooper pair reside in different electrodes. I propose here a situation in which 
the spin-cr electron making the Cooper pair is forced to couple to a vector potential (by propagating 
on a loop in the presence of a magnetic flux) while the spin-(— a) electron is not directly coupled 
to the magnetic flux (it propagates in an ordinary ferromagnet without any hole). As I show, the 
spin-(-a) current oscillates as a function of the flux coupled to spin-a electron. This new type of 
Aharonov-Bohm effect can be viewed as a direct consequence of non separable correlations, and may 
be tested in future experiments. 

The article is organized as follows. As a warm up exercise, I first calculate the transport formula 
associated to non local Andreev reflections in section ^. I introduce the effective single site formalism 
in section ||, and provide a discussion of higher order Feynman diagrams. Transport of linear super- 
positions of correlated pairs of electrons is derived in section |||. The Aharonov-Bohm experiment is 
discussed in section [|. I present in section ^ a simple model for superconducting propagation along 
a domain wall, which appears to be a problem closely related to non separable correlations. Final 
remarks are given in section |7[ 

2 Warming up: transport of non local Cooper pairs 

I start by considering the model presented on Fig. |l] in which a superconductor is connected to N 
ferromagnetic electrodes. I first determine the superdoncuting gap in a self-consistent way, and next 
calculate the current I n flowing into one of the ferromagnetic electrodes. 

2.1 Effective single site Green's functions 

Let us first describe the superconductor and ferromagnetic electrodes on the basis of effective single 
site Green's functions. The Nambu representation of the superconductor Green's function takes the 
form [21] gx^iuj) = fl + gcr x , with 



g(u ± in) = Tip N —= u;± ' l ' n ; an d ± ^) = 7rpN . ^ 



^A 2 - (lo ± in) 2 ' ~ ' ^A 2 - (u> ± in) 2 

where p^ is the normal density of states. The Nambu representation of the effective single site 
ferromagnetic Green's function reads ga' k R ak = ±iw \pk,l[I + & z ]/2 + Pk,i[I — <5' z ]/2|, with pk >a the 
spin-a density of states in the ferromagnetic electrode ctk ■ The Nambu representation of the hopping 
matrix element is t x ,u h = tx,a k & z - Using effective single site Green's functions means that the 
superconductor is viewed as as being "zero dimensional" . This is a valid assumption if the distance 
between the contacts is smaller than the coherence length. This is precisely the condition under which 
non local Andreev reflections occur [jl2|J. Therefore, representing the superconductor by effective 
single site Green's functions is a self consistent assumption. 
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Figure 1: Schematic representation of the model considered in section || in which a superconductor 
(site x) is connected to N ferromagnetic electrodes (sites a±, ... , ajy)- A voltage V is applied on the 
superconductor. 



2.2 Self consistent determination of the superconducting gap 

To determine the superconducting gap, we need to calculate the Gorkov function. First, we solve the 
Dyson equation 



G 



R.A 



N 



9x'x ^2 



x,a k 9a k ,a k L ce k ,x 



R.A 



tr 



k=l 



q R > A 
Hx,x • 



(1) 



The relevant parameters appear to be the spectral line-width associated to spin-<7 electrons: T a = 
J2k=i r k,a, with T kj(T = (t ak ,x) 2 Pk,a- Solving Eq. pleads to 

G* x = ^{[g + i7r(f 2 -g 2 )]i + fa*}, (2) 

with T> = 1 — iirg(Ff +7r 2 (/ 2 — g 2 )T^Ti. To calculate the superconducting order parameter, we 
need to solve the Dyson-Keldysh equation G H = (I + G R (g> W) ® g + ~ ® (I + W ® G A ), where the 
convolution includes a sum over the labels x and a^. Noting X a = (1 — ingT-^/T), Y a = infT^/T), 
and using Eq. |2[ we obtain the exact expression of the Nambu component of the Keldysh Green's 
function: 



[G. 



x,x\ 2 A 



= 2iimp{uj) x 

{ Pg (x T F T + y l x{) + Pf {x{x i + y ; y T ) 
■ 1 y 1 (x^-i\ + -L,yUx 1 - 



(3) 
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where np{oj) is the Fermi distribution, and I used the notation p = p g I + pfd x = Im[<7 A ]/7r. The 
superconducting gap is obtained via the self-consistency equation 21]: 



p+oo 



A = U 



/ duj/{2m)[G + -{u)] 2 . l 



J — oo 



with U the microscopic attractive interaction. The dominant contribution arises from the large-|u 
behavior, which leads to the BCS-type relation: 



with D the bandwidth of the superconductor. As an example, I consider a coupling to two ferromag- 
nets. With a parallel alignment of the magnetizations in the electrodes, one has I\ = 27 and r\ = 0. 
With an antiparallel alignment, one has Tt = T \ = 7. The ratio of the two gaps is found to be 



which shows that the spin polarized environment generates a reduction of the superconducting gap, 
stronger with an antiferromagnetic alignment than with a ferromagnetic alignment. As a consequence, 
the transition temperature of the superconductor is larger if the electrodes are in a parallel alignment. 
This is valid whatever the transparency of the contacts. 

To make contact with possible experiments, one should keep in mind the following points: 

(i) Eq. H should be strictly speaking applied only if the dimension of the superconductor is small 
compared to the coherence length. 

(ii) It is nevertheless possible to extrapolate to the physics occurring when a bulk superconductor 
is connected to ferromagnetic electrodes. In this case there is a non homogeneous gap in 
the superconductor, which is reduced at the interfaces with ferromagnetic electrodes. If two 
ferromagnetic electrodes are at a distance smaller than the coherence length, the local gap 
should follow qualitatively Eq. ||. A transport experiment involving two electrodes at a distance 
smaller than the coherence length would be sensitive to the local gap. 

In the remaining of the article, I focus on transport properties. While keeping in mind that the 
superconducting gap should depend on the spin orientation of the ferromagnetic electrodes, I will 
not determine systematically the superconducting gap in a self consistent manner. Instead, I assume 
that the gap is a fixed quantity, but I could easily inject a posteriori the self consistent value of the 
gap in the transport equations. 

2.3 Non local Andreev reflections 

I use Green's function techniques to evaluate the current I n flowing into the ferromagnetic electrode 
OLiyi CIS cl trace over Nambu space: 



Since it is trivial to restore the quantum of conductance in the transport formula, I assume that 
e/h = 1 throughout the article. 




(4) 




(5) 
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The Hamiltonian decomposes into two contributions: 



N 



k=l 



where 7i x and H ak are associated to the superconducting and ferromagnetic electrodes, and W is 
the coupling term: 



TV 



VV — ^afc,^ [ c a k c x + c x c a k ] 



fc=l 



It is a standard method to treat W as if it were a perturbation, and to sum to all orders the resulting 



perturbation theory. The final transport formula is then non perturbative in W [22|. Therefore the 
method is well suited to describe the physics of high transparency contacts. 
The chain of Dyson-Keldysh equations associated to Fig. |l| takes the form 



N 



N 



A 



k=l I 1=1 

tx,a„ \l + G an X t x>an 

^ ' tx,a n Ga n ,xtx,a k 9a k ,a k ta k ,xG 
k^n 



?<+- t r A 

yOtntCHn <*n,X ^*X,X 



A 

x.x • 



(7) 

(8) 
(9) 



The spin-cr Andreev current flowing into electrode a n is found to be 

= 47r 2 r„ j(T r_ (T J du [n F (uj - eV) - n F (uj)] \G£ )X>1j2 \ 



(10) 



The Andreev current Eq. 10 is similar to the one of conventional Andreev reflection [21|, except for 
the density of state prefactors appearing in the spectral line-widths. These prefactors reflect the fact 
that a spin-<7 electron can be transfered in electrode a n only if it paired with a spin-(— a) electron 
either in the same electrode (in which case the Cooper pair is a local object), or in another electrode 
(in which case a non local Cooper pair is transfered). This shows that the Andreev current is built 
up from a sum of all possible non local Cooper pair contributions, and generalizes the situation 
considered in Ref. JO]]. 



3 Effective single site Green's functions formalism 

Now I describe the method that will be used throughout the remaining of the article. The strategy is 
to use Green's functions techniques in an effective single site formalism. This method has proved to 



be extremely powerful in the context of superconducting quantum point contacts [21|. The advantage 
is that simple algebraic expressions can be manipulated without using spectral representations. The 
basic assumption is to consider that a given site represents a phase coherent continuum of states. I 
first take the example of the metal - metal - metal junction and show how to interpret the resulting 
transport theory. Since I consider in this section circuits without superconducting elements, the 
Nambu space structure is diagonal, and I can use spinless fermions without any loss of generality. 
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Figure 2: Notations for the Green's functions. In (a): no link is disconnected; the corresponding 
Green's functions are noted Gj.j. In (b): the link a - a is disconnected; the corresponding Green's 
functions are noted gij. In (c): The two links a - a and b - (3 are disconnected. The Green's 
functions are noted hij. One has g a>a = h ata . 

3.1 Treating all couplings on the same footing 

Let us start to calculate the transport formula corresponding to the metal - metal - metal junction 
shown on Fig. [2| We note G the Green's function of the fully connected system, g the Green's 
function with t a ^ a = 0, and h the Green's function with t a , a = tb,p = (see Fig. |2|). 

I treat simultaneously the two couplings t ay0l and i& « on the same footing. Namely I go directly 
from the G's to the h's without ever using the g's. The chain of Dyson equations for the Green's 
function takes the form 

G a ,a — ha,a ~i~ h a ^ a t a ^ a h a ^ a t a ^ a G a a + h aa t aa h a ptpbGb^a (H) 
Gb,a = hb,btb,php % ptp t bGb ta + hb^tb,php %a ta,aG a ,a- (12) 

The effective single site Green's functions are ha,' R = ±i7rp a , hf'- R = ±iirp (i, j = a,/3), hf' b R = ±iirpb. 



I use the notation 7 a = Tr 2 \t a , a \ 2 pa and 7fe = ^ 2 \tb,p\ 2 Pb- Eqs. 11 - |i~2] are solved in a straightforward 
fashion. For instance, one obtains 

l + P(la + lb) 

Next, I calculate the Keldysh component 

Ga,°> = C''" G a ,aTa,a\ h^ a t a>a G a a + G a a t a ^ a h^ a [l + ta,aG a Q \ + G a a t a ^ a h a ptp^bG\ ) OL (13) 
+ G a,btb,l3hp^ [l + t«,aGa, a ] + G^btb^hp^tp^G^ + G^ptpfih^t^pGp^. (14) 

It is straightforward to evaluate the six terms to obtain the transport formula 

l a ,a= I au- ; --o - eV) - n F (io)\ . (15) 

J [1 + P{7a + 76) 
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Figure 3: Representation of the current flow. Iq is a disconnected current and only Ip b is relevant. 



3.2 Decimating part of the circuit 



Before discussing the physics of Eq. 15, I give another derivation of this transport formula. The idea 
is to use first a perturbation in t b a and, in a second step, make a perturbation in t aa . Namely, part 
of the circuit is decimated exactly, and I look for an effective theory for the remaining sites. 
The perturbation in t a ^ a leads to 

^a,a = ^a,a^a,a9a,a [l + ^ot,aG a ,a\ + \\ + ^a,a^a,a] dt,a ^a,aG aa - (16) 

Therefore, I need to determine g^~ and g^~ in terms of the h's. The chain of Dyson equations reads 



9a,a = h at a + h a ^tpfi9b,a 
9b,a = hb,btb,f39l3,a 
9/3, a = h/3 t a + hpfitpjg^a- 



(17) 
(18) 
(19) 



The equation for the Keldysh component is 

= 9a,b t b,/3h'pfitp tb g ba + h+~ + h^t/3 ib g ba + g^bh^h^' + ga t /3t(3,bh b btb,/3gp, a - (20) 
It is straightforward to solve Eqs. |l7| - 2C, inject the propagators into Eq. [l6| and find the current 



Eq. 15. In this approach, the effect of the coupling 0— b is just to renormalize the density of states 



at site a: g A 



2mn F {uS)p a , g 



A,R 



±i7rp a , with 



Po 



1 + Plb ' 



(21) 



where p is the density of states in the intermediate region and j b = 7r 2 |t{,,/3 | 2 j°6- We will apply this 
decimation method latter in more complicated situations. 

3.3 Interpretation of the transport formula 

The current Eq. |15] can be decomposed into two contributions: I a ^ a = Iq + Ip >b , with 



duo- 



[1 + p{la + lb)\ 



[uf(uj — eV) — uf(uj)] 



(22) 
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Figure 4: The diagrams of order p 2 ^ a lb contributing to Ipj,. The top diagram has a unit charge 
(Q = 1). The bottom diagram has a zero charge (Q = 0). 



and 



du- 



[hf(lo — eV) — np{oj)] . 



(23) 



[1 + P(la + lb)Y 

Is is straightforward to calculate directly the current flowing through the link /3 - b, and show that 
it coincides with Eq. |23|. Therefore, we are lead to consider that the region a - (3 can absorb part 
of the electrical current (see Fig. |3|). Only is the transmitted current Jg^ a quantity of interest. To 
understand the appearance of the fake term Iq, it is useful to notice that an effective single site plays 
two roles at the same time: 

(i) It transmits current to the rest of the circuit. This is why there is a finite contribution Ig^. 

(ii) It represents a continuum of energy levels; therefore, it also plays by itself the role of a reservoir. 
This is why there is a finite current Iq. 



Now the expression (23) of Ipj, contains also a renormalization prefactor which generates an 
infinite series of Feynman diagrams. It turns out that disconnected contributions are also present in 
these diagrams. To obtain the final form of the physical current, we expand Eq. |23| in an infinite series 
of Feynman diagrams, and operate on this series with an operator Q that multiplies each Feynman 
diagram its charge: 



'phys 



dujQ 



lalb 



[^+p{la+lb)Y 



[uf{u> — eV) — np(w)] 



(24) 
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Figure 5: The four diagram of order p^ialb contributing to I$h- The two top diagrams has a unit 
charge (Q = 1). The two bottom diagrams have a zero charge (Q = 0). 



In the case of the metal -metal -metal junction, Eq. |24| should coincide with the transport formula 
obtained by Caroli et al. 



'phys 



/dui — 
1 



P 2 lalb 



+ p(la+lb)\ 



[uf(^ — eV) — np(uj)] 



(25) 



To verify this, I calculate the lowest order diagrams associated to Eq. 24. This will give us the 
opportunity to show how the operator Q operates. 

3.4 Feynman diagrams 



To count how many diagrams appear at a given order, we expand Eq. 23 in a power series: 
2p 2 lalb 



(l + p(7a+7&)) 5 



2/0 7a7fe - 4 P 7a lb (la + lb) + €>P lalb (la + lb) + 



(26) 



From what we deduce that there are two diagrams at order 7 a 7b, four diagrams at order 7„7b, six 
diagrams at order 7„7&, twelve diagrams at order l\l1-, etc ... Each of the diagrams corresponds to a 
possible contraction of the density of states. For instance the two diagrams at order 7 a 7b are shown 
on Fig. ||. The four diagrams at order 7„7b are shown on Fig. [5|. The twelve diagrams at order 7^7^ 
are shown on Fig. [||. Using these diagrams to evaluate Eq. |24|, we obtain the expansion 



Vys = J ^ [P 2 lalb ~ 2p 3 lalb(la + 7b) + ^P^lalbda + Ibf + ■■■} X [n F (uJ - eV) - n F (u)} 

which coincides with the expansion of Eq. p^. 



9 



00 oo 



o^ 

o^ oo 



oo^o 

00 ^0 








00 00 







0^00 



00^0 




Q=l 



0= 
0= 


















C3> 



Q=2 ; 



> Q=° ; 



Figure 6: The twelve diagram of order p Jalb contributing to Ipjy. The charges are indicated. 

4 Transport of linear superposition of correlated pairs 
4.1 EPR paradox 



It is worth recalling the well known EPR paradox [23]. Let us consider a collision of two spin- 
1/2 electrons. The collision is represented by an antiferromagnetic exchange interaction J(r)Si.S2, 
where J(t) is constant in the time interval [0, t] and zero otherwise. After the collision, the system 
is represented by the EPR state 



cos (Qt) \ + -) -isin(Oi)| -+), 



(27) 



with fi = J/2. Let us imagine that two observers 0\ and O2 use an apparatus (i.e. a Stern and 
Gerlach magnet) to measure the spin orientation of Si and S2. The observator 0\ measuring Sf 
would find 1/2 with a probability cos 2 (Qt) and —1/2 with a probability sin 2 (fit). Once the observer 
0\ has made his measurement, the wave function collapses and with a probability one, the observer 
O2 measuring S| would find Sf = 1/2 if Oi has measured Sf = —1/2, and Sf = —1/2 if 0\ has 
measured Sf = 1/2. Quantum mechanics contains non separable correlations, which is the famous 



EPR paradox |2^]. As pointed out by Bell, there exists a basic difference between quantum mechanics 
and hidden variable theories Q . It has been demonstrated experimentally that Bell inequalities were 
violated with photons ||, therefore ruling out hidden variable theories (at least with photons). As 
already mentioned in the Introduction, it is an open question to test the EPR paradox with electrons. 
We do not address this issue directly in this article. Nevertheless, we find an interesting connection 
between transport properties and quantum measurement (i.e. projections) of linear superpositions 
of correlated pairs. 
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Figure 7: Schematic representation of the model considered in section ||. The site x is supercon- 
ducting. The intermediate regions a - a' and /3 - (3' are ferromagnetic. The electrodes a and b are 
ferromagnetic. 



4.2 Fabrication of linear superpositions of correlated pairs 

Now I use the formalism presented in section || to discuss the physics of linear superpositions of 
correlated pairs of electrons. I consider that a ballistic ferromagnetic region otk - a' k is inserted in 
between the superconducting site x and the external ferromagnetic electrodes (see Figs. [?] and |||). 
To simplify the discussion, I make the following assumptions: 

(i) There are no pair correlations between the bottom electrode on Fig. |7| (used to impose the 
superconductor chemical potential) and the other ferromagnetic electrodes. This means that 
the corresponding contact on the superconductor is separated from the other contacts by a 
distance much larger than the superconducting coherence length. 

(ii) All ferromagnets are fully polarized. Therefore, there is no transmission of local Cooper pairs. 

(iii) The ferromagnetic intermediate regions are phase coherent. Present time technology does not 
allow the fabrication of such devices. Nevertheless, the model is interesting per se. 

Our analysis of the situations presented on Figs. and ||] is the following. The wave function 
in the ballistic regions is a linear superposition of correlated pairs. For instance, in the presence of 
three fully polarized ferromagnetic regions connected to the superconductor having a spin orientation 
a a =|, op = (Tj =[ (see Fig. ||), the wave function is 



In the presence of spin orientations a a = ap = j, cr 7 



+ 



y t/3P/3,l + t 1 p lA 
I, we have 



|0>. 



(28) 



+ 



taPa,] 



|0>. 



(29) 



These wave functions are the only ones that can guarantee the correct pairing associated to the 
formation of correlated pairs of electrons. The coefficients of the wave function can be obtained by 
evaluating the Gorkov function. With the configuration a a = op =|, cr 7 =|, we find 



1,2 



7T t a ^tyP a (J3)^P^l 



1,2 



From what we deduce Eq. 28. We will come back in section 4.4 on the meaning of this wave function. 



The dictionary between the Cooper pair wave function Eq. 
Eq. ^ is the following: 



+ 



+ 



and the spin-1/2 wave function 

(30) 
(31) 
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Figure 8: Representation of a device in which three ferromagnetic electrodes a, (3, 7 are connected 
to a superconductor. The additional electrodes a, b, c are used to to perform a measurement of the 
linear superposition of correlated pairs. 

The relevant question that should be asked now is to determine whether it is possible to make 
a quantum measurement of the linear superposition of correlated pairs, similar to the Stern and 
Gerlach measurement of the spin-1/2 EPR state. Namely, we want to interpret transport theory in 
terms of projection operators. 



4.3 Transport formula 

The strategy to solve the model on Fig. [?] within the effective single site Green's function approach 



is to use the two step perturbation theory presented in section 3.2. I first treat the link a 1 - a 



in perturbation (see Fig. 0). Since the superconducting site has been disconnected and there is 
no source of spin flip scattering, the Dyson equations decouple into a spin-up and a spin-down 
component. We note hij the Green's functions of the disconnected system, i.e. with t a ^ k = 0. The 
Dyson equation for the Keldysh component has already been given in section ||. The solution takes 
the form g+~ a = 2innF(^)pa,a, and ga'a,a = H^Pa,cr , with the renormalized density of states 



Pa,a 



Po 



1 + 7T |ta, a| Pa,aPa,a 



Next, I use a perturbation in t x ,a k to calculate the Andreev current. This is done by noticing that 
the physical current I a<a is a finite fraction of the current incoming at site a (see Eqs. 



2D: 



Fa,Q!| PaPa 
1 + ^ 2 \ta,a\ 2 PaPo 



(32) 
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Figure 9: The relevant lowest order Andreev reflection diagrams, (a) is of order Pa n Pa m {Pa n ) 2 {Pa m ) 2 ■ 
(b) is of order p an (pa m ) 3 {paj 2 (pa m ) 4 ■ (c) is of order {p a „f pa m {PaJ 4 (pa m ) 2 ■ The dashed line 
represents the anomalous propagation in the superconductor. 



I deduce from Eq. 32 that the Andreev current Eq. [T(] is the sum of all possible Cooper pair trans- 



missions: 



where 



J n,a = 471-2 J duJ [ n H w - eV) - n F (u)) ^ Q [u n ^ a 'U m \G^ x ,i. 



21 ) 



— 2 lf 1 2 r - |2 



Pa n ,aPa n ,(J 



+ ^ 2 K n ,a' n \ 2 Pa n ,aPa n ,aY 



(33) 

(34) 



The action of Q in Eq. 33 has already been given. The form Eq. 33 of the current can also be obtained 
by making directly a perturbation in t x ^ ak to evaluate the current through the link a n — a n . The 
density of state prefactors appearing in the lowest contribution in Eq. |33| corresponds to the diagram 
on Fig. H (a). 



4.4 Quantum measurement of the superposition of correlated pairs 

Let us now interpret the transport formula Eq. |3^ in terms of quantum measurement of the superpo- 
sition of correlated pairs. We assume that the external ferromagnetic electrodes are weakly coupled 
to the intermediate regions. As a consequence, we make the approximation p a = p a . The interfaces 
with the superconductor can be either tunnel or high transparency contacts. The general form of 
the linear superposition is 

W) = N- 1 ' 2 VW^,t/^,l<,T<J°>' (35) 

P,9 
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where 

p q 

It is useful to define a projection operator associated to the correlated pairs of electrons (p f, q [): 

The spin-up current through electrode a p can be rewritten in the form 

/ p = 4vr 2 ^ 2 ^ I ^[^(a;-eF)-n F (a;)]7^ T 7 09ji |^|P M |V')|V^,i 1 2| 2 ) (36) 

q J 

where 7a p{9) ,t = ^ rapfa) I /%(<,), t(i) denotes the spectral line-width associated to the external elec- 
trodes. We first notice that the current is a sum over all possible Cooper pair transmissions, which 
we already found in section ^. What is new here is that there is one intermediate quantum state is 
involved (the linear superposition). The transport formula can be expressed as suitable projections of 
this quantum state. These projections correspond to the physical processes by which current is car- 
ried away from the quantum state, namely the transport of correlated pairs of electrons. Therefore, 
the whole picture is consistent: we first identified a wave function, and next showed that transport 
can be interpreted in terms of quantum measurement of the wave function. In our opinion, the 
form ( |36| ) of the transport formula is not specific to our model, but is valid independently on the 
microscopic origin of the pairing potential. 



This behavior is also consistent with general considerations on decoherence [^4|, 25, 26, p7[ |, at 
least on a heuristic basis. I conjecture that there are three distinct stages involved in the decoherence 
mechanism, which occur on different time scales: 

(A) The initial quantum superposition Eq. being a Schrodinger cat of correlated pairs of elec- 
trons. 

(B) An intermediate stage in which the Schrodinger cat superposition is lost but there are still pair 
correlations. Namely, the spin-up (-down) electrons have tunneled into the specific reservoirs 
p (q). The correlated pair still behaves like a two-electron phase coherent object. 

(C) The final stage in which the two electrons making the correlated pair behave like independent 
quasi-particles, with no more pair correlations. 

The transition (B) — > (C) is usually known as the proximity effect. I conjecture that the projector 
appearing in Eq. |36| is the signature of the decoherence process (A) — > (B). Note that because two- 
electron correlations are preserved until stage (C), there is an interesting physics g omg on at stage 
(B) which is explored in the next section. 

The reader might object that the intermediate regions may contain thousands of electrons. It may 
then appear counter-intuitive to propose Schrodinger cat made of superpositions of a large number 
of electrons. This is certainly true, and one should keep in mind that the wave function Eq. is 
by no way the many body ground state wave function. Instead it is intended to describe only the 
degrees of freedom that participate to transport. 
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Figure 10: Geometry of the Aharonov-Bohm experiment. The current I flowing into the left electrode 
is modulated by the flux enclosed by the other electron making the correlated pair. The non separable 
correlations between the left electrode and the Aharonov-Bohm loop are represented schematically. 



5 Aharonov-Bohm effect 

Now I propose a device that can be used to probe non separable correlations. For this purpose, 



I consider the geometry on Fig. 10, in which one side of the superconductor is connected to an 
Aharonov-Bohm loop (being made of a normal metal), while the other side is connected to a normal 
metal that does not contain any hole. I show that the current flowing through the electrode with no 
hole is modulated by the flux enclosed by the normal metal loop. This means that the current due to 
one electron making a correlated pair is modulated by the flux enclosed by the other electron, which 
is a manifestation of non separable correlations. The electrodes are considered to be normal metals 
in which phase coherence can propagate over large distances. 



I consider the circuit on Fig. 11 that is used to model the situation on Fig. |10[ The hopping 
parameters are complex number, for instance i a>a / = t exp (i(p a ,a') = texp[2i7T(f>/(4(j)o)], and t a * A = 
texp (if a ', a) = texjp [— 2i7r^>/(40o)]. Since I use effective single site Green's functions in which a' and 
a" are represented by the same single site, this choice of the hopping insures that § (fij = 2ir(j)/(f>Q. 
The first task is to solve the chain of Dyson equations 

G a ,a = 9a,a + 9a,ata,a'G a i t a + g a) ata,a"G a " t a (37) 

G a ',a = ga',a'ta',aG a ,a + 9a',b'tb',bGb,a (38) 

G a ",a = g a ",a"ta",aG a ,a + g a ",b"h",bGb,a (39) 

Gb,a = 9b,bh,b'Gb' ,o + 9b,bh,b''Gb" ,a (40) 

Gb',a = 9b' ,b'tb' ,bGb,a + 9b' ,a't a ' , a G ' a ,a (41) 

Gb",a = 9b",b"tb",bGb,a + 9b",a"ta",aG a ,a- (42) 

Next we need to calculate the Keldysh component 

G t t a = [ 1 + G a,a t a',a + G^ >a nt a >f ia ] [l + i a ,a'G^a + W" G a",a] ( 43 ) 

+ [ G ^6"V,& + Ga y b'h',b\ 9b.b [tbfi'Gy a + hpiGy, a ] (44) 
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Figure 11: The circuit used to represent the situation on Fig. 10. x is superconducting, (3 is ferro- 
magnetic and the remaining sites are normal metals. 



+ 


" a,a L a,a' 9 a i a i ^a' ,a<J a,a 


(45) 


+ 


(~<R i H — j. /~iA 
^ a,a l a,,a" il a" ,a" la " , a(Jr a,a 


(46) 


+ 




(47) 


+ 


Ga,btb,b' 9 t', a ' *<*' ,a^a,a 


(48) 


+ 


Ga,ata,a" 9a", b" ,b^b,a 


(49) 


+ 


G a,b t b,b"g^a" t a",a. G a,a, 


(50) 


+ 


Ga,btb,b' 9b< ,b' tb ' ,bGb t a 


(51) 


+ 


Ga,b t b,b"9b'^b" t b",bG b A ',a- 


(52) 



One finds Ga'^ = ±z7rp a , and G+ a = 2mn F {uj)p a , with the renormalized density of states 



Pa 



p a [l + p^y + y')] 



1 + (Pa + Pb)(l' + 7") + ^I'i'paPb [1 " COS (27T0/0Q)] ' 



and where I used the notation 7' = ir t p' and 7" = n i p" with f = |t a 



1*6,. 



(53) 

1*6,6" I • 



Is is straightforward to use the same procedure as in section |4.3| to eliminate disconnected pro- 
se- 
form 



cesses, which leads to if = j£ local + i£ Tinanlocal + ^ iiIlon i oca i, where the local contribution takes the 



^locai = 4vr \tx,p\ PMPM J du[n F (u-eV) - n F (u)] \G~ X ^ 2 \ , (54) 
and the non local contribution is 

^ 2 %A 2 %A 2 p^p\- a Pb,-a [(7-a + i-af + ^'-vl-M ~ COS (2vr0/0 O )]] ' 



/3,<r,nonlocal 



duiQ 



[1 + ( Pa ^ a + P b,-a){i- a + 7- J + 2Y-al-aPa,-*Pb,-v[l " COS (2lT <f) / (j) )}] ' 



[n F (uj - eV) - n F (uj)} \G X - X 



21 > 



(55) 



where the action of Q has already been given. The non local current Eq. |55j is periodic in 4>/4>o. Such 
oscillations are a signature of non separable correlations. 

The origin of these oscillations can be understood in simple terms in the limit i a ,a;> *6^ <C 1, where 
the physics is the one of resonant tunneling. There are discrete energy levels on the loop, the energy 
of which depends on the magnetic flux. Single electron transmission is modulated by the flux because 
the transmission is larger when there is one resonant level. Now the Andreev current is carried by 
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electrode 



( no resonant level ) 




Superconductor 



Figure 12: The diagram involved in the absence of resonant tunneling in the right electrode. This 
diagram does not carry electrical charge. 



correlated pairs of electrons. If the Aharonov-Bohm loop is off-resonance, one electron making the 
correlated pair (for instance a spin-up electron) cannot be transmitted. As a consequence, the whole 
correlated pair cannot be transmitted. There is therefore no spin-down current. 

This situation is not far from the behavior of the multichannel ferromagnet - superconductor 
junction. It is well known that in such junctions there is no current transmitted in the channels 
having only a spin- up Fermi surface [p8fl . 

The following equivalent explanation can be found. There is one virtual state in which one Cooper 
pair has been extracted from the superconductor, one electron is in the left electrode and the other 
electron is in the right electrode. Since there is no resonant level, the electron in the right electrode 
is backscattered onto the superconductor, undergoes an Andreev reflection in which one Cooper pair 
is transferred into the superconductor, and a hole is transfered in the left electrode. The whole 



processes is represented on Fig. |12j. A quick look at this diagram shows that there is no current 
transmitted. Note that the processes involved in the multichannel superconductor - ferromagnet 



junction can be interpreted with the same diagram [29|. 

We stress that these two pictures based on resonant tunneling hold only if the interfaces have 
low transparency contacts. Nevertheless, since we found a systematic way to handle higher order 
Feynman diagram, the model is valid also for high transparency contacts, in which case Aharonov- 
Bohm oscillations are also predicted. 



6 Propagation of superconducting correlations along a domain wall 

Now, I illustrate the consequence of the model regarding the propagation of superconducting cor- 
relations along domain walls in ferromagnets. Such propagation can generate an enhancement of 
the proximity effect in ferromagnet / superconductor heterostructures, not against recent experi- 
ments [13, 14, [l5|, 30 1 . We use the geometry on Fig. |l^ to address this problem in an effective single 
site formalism. If site a has a spin-up orientation and site has a spin-down orientation, there is one 
domain wall in the junction. In the presence of partially polarized ferromagnets, a spin-u electron 
can be transferred across the domain wall because of the hopping matrix element t'. This provides 
a minimal model for Cooper pair penetration along domain walls. This model in which the domain 
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a 




P 

Figure 13: The circuit used to represent the propagation of superconducting correlations along a 
domain wall. There is a domain wall if a has a spin-up magnetization and (3 has a spin-down 
magnetization. I note t = t a , a = ta,p = %,a = tb,p and t' = t a ^. Site x is superconducting. 



wall has a vanishing width is expected to be relevant for describing a situation in which the width of 
the domain wall is smaller than the superconducting coherence length. 



The Dyson equations corresponding to Fig. 13 take the form 

G a ,a — 9a,a ~i~ 9a,ata,aG a ,a ~\~ 9a,ata,fiG (3,a 
G a ,a = 9a,ata,aG a ,a + 9a,otta,pG f3,a + 9a,aTa,b 

G h 

G/3,a = 9p,^t/3,aG a ,a + 9l3,j3tj3,aG a ,a + 9l3,!3tf3,bGb,a 
Gb,a = 9b,btb,aG a ,a + gb,btb,f3Gf3, a , 



(56) 
(57) 
(58) 
(59) 



A,R 



and similar equations hold for the Keldysh component. These equations are solved into G a ,a 
±z7T/5 a , G^~ = 2iirnF(ui)p a , with the renormalized density of states 



Pa = Pa 



1 + TT 2 t' 2 p a pf3 + TT 2 t 2 p b [p a +Pf3 + 2mt' Pa pp\ 
1 + ir 2 t' 2 p a P/3 + n 2 t 2 (p a + Pb) [Pa + P/3 + 2int' p a pp] ' 



(60) 



I deduce from Eq. 6C that the Andreev reflection current takes the form 



4vr 6 i 4 t 8 



duiQ 



n.n 



Px,aP a ,aPb,c 



(pa,a + PP,oY + ±-K 2 t' 2 p 2 ^y p 



[uf(oj — eV) 



+ ir 2 t' 2 p a:(7 pi3,a + vr 2 t 2 [p a!(T + pp t(T + 2i7Tt'p ata p f3:(7 } 



x (61) 



(62) 



In the limit of fully polarized ferromagnets, the current Eqs. 61, ^32] vanishes if the sites a and /3 
have the same spin orientation. In the presence of a domain wall, a and (3 have an opposite spin 
orientation, in which case Cooper pairs can be transfered across the junction even in the presence of 
fully polarized ferromagnets. 

The effect proposed here is related to two other well known problems: 

(i) The it -junction: It is well known that Cooper pairs cannot propagate over large distances in a 
ferromagnetic metal. In fact, they cannot even propagate at all in a fully polarized single domain 



18 



ferromagnet. In partially polarized ferromagnets, Cooper pair penetration gives rise to the 7r- 



j unction physics that has been investigated by theorists since several decades [31, |33|, gj], an d 



has been obtained experimentally only recently [35|, 36]. We have shown here that Cooper pairs 
can propagate in fully spin polarized ferromagnets if there are domain walls in the ferromagnet. 
The number of conduction channels is equal to the number of domain walls. 

(ii) Cryptomagnetism: It has been established a long time ago by Anderson and Suhl that fer- 
romagnetism and superconductivity can accommodate each other in the same system if the 



ferromagnet acquires a domain structure |37|, |3g]. Our proposal can be viewed as the proximity 
effect version of cryptomagnetism. 

7 Conclusion 

To conclude, I have provided a minimal formalism to treat non separable correlations generated when 
a superconductor is in contact with ferromagnetic or normal metal electrodes. It has been shown 
how to solve the problem of disconnected contributions arising in this formalism, and it has been 
given a systematic way to handle the infinite series of Feynman diagram. 

We answered the four questions given in the introductory section, namely: 

(i) The superconducting gap depends on the relative spin orientation of ferromagnetic electrodes 
connected to a superconductor. This might be probed in experiments. 

(ii) It is possible to fabricate linear superpositions of correlated pairs of electrons. Transport in 
these systems can be expressed in terms of projections of the linear superposition of correlated 
pairs of electrons. The wave function is projected on the relevant quantities that participate to 
transport, i.e. the correlated pairs of electrons. 

(iii) I proposed a new Aharonov-Bohm experiment to test non separable correlations. In this exper- 
iment, the spin-ex electron making the Cooper is forced to couple to a magnetic flux while the 
spin-(— a) electron has an ordinary propagation. It has been shown that the spin-cr current 
oscillates as a function of the magnetic flux coupled to the spin-down electron. It has been 
predicted that the effect appears with low and high transparency contacts, and independently 
of the spin polarization of the electrodes. Therefore, normal metal electrodes can be used, in 
which phase coherence can propagate over large distances. 

(iv) It has been shown that superconducting correlations can propagate along ferromagnetic domain 
walls. The model is strictly speaking valid if the width of the domain wall is small but we 
are confident about forthcoming generalizations of the model. The number of Cooper pair 
conduction channels is equal to the number of domain walls. This may apply to existing 



experiments in heterostructures [13, 14, 15 1 which have not received a satisfactory explanation 



up to now in spite of various theoretical attempts [15, 30]. 



To end-up let us mention that the work presented in this article calls both for new experiments 
and for more sophisticated theoretical descriptions. It is in general needed to incorporate realistic 
ingredients in the models, which will be the subject of future publications. Nevertheless, we have no 
doubt that the same physics will be obtained in the presence of realistic constraints. 

The author acknowledges fruitful discussions with P. Degiovanni, D. Feinberg and M. Giroud. 
This work is supported by the French Ministry of Research under contract ACI 2086 CDR2. 
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